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ABSTRACT 


thermal  conductivities  of  Praseodym^^ijm-Cer  ium-Oxide 
(PCO)  and  Saran  Carbon  have  been  experimentally  investigated 


using  a  steady-state  heat  transfer  technique.  The 
investigated  substances  are  used  as  adsorbents  in  adsorption 
compressors  being  developed  for  spaceborne  refrigeration 
applications.  The  objectives  of  the  investigation  were  to 
determine  the  thermal  conductivities  and  establish  their 


temperature  dependency.  Data  were  collected  for  the  PCO 

v .  j/s  Jy  y 

over  a  temperature  range  of  300*0!  to  600*  C ,  and  /)*  C  to  200  C 

for  the  Saran  Carbon.  The  thermal  conductivities  were  found 

to  have  a  strong  temperature  dependency.  In  particular,  the 

results  for  the  PCO  showed  a  temperature  dependency 


indicative  of  some  thermal  radiation  effects. 
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CHAPTER  1 


INTRODUCTION 

In  space,  a  refrigeration  system  is  needed  to  provide 
cooling  to  many  instruments  located  onboard  the  spacecraft. 
Typically,  refrigeration  systems  have  relied  on  mechanically 
driven  compressors  to  compress  the  working  fluid.  These 
compressors  are  not  advantageous  for  space  application 
because  they  have  wear-related  moving  parts  which  limit  the 
mission  duration.  An  alternative  is  to  replace  the 
mechanically  driven  compressors  with  adsorption  compressors. 
Because  of  the  lack  of  moving  parts  non-mechanical 
adsorption  compressors  offer  a  longer  operating  life,  far 
exceeding  mechanically  driven  compressors.1 

Adsorption  compressors  are  essentially  canisters  filled 
with  a  porous  material  that  desorbs  and  adsorbs  gases  as  the 
temperature  in  the  canister  is  increased  and  decreased.2  By 
thermally  cycling  a  number  of  canisters  filled  with  a  porous 
adsorbent,  it  is  possible  to  provide  a  flow  of  high  pressure 
gas  for  the  refrigeration  system. 

Currently,  the  Jet  Propulsion  Laboratory  (JPL)  in 
Pasadena,  California  is  developing  a  laboratory  model  of  an 
absorption  cryocooler  system1’3  which  uses  Saran  Carbon,  and 
Praseodymium-Cerium-Oxide  (PCO)  for  the  porous  adsorbent. 

In  conduction  with  JPL,  numerical  studies  are  being 


conducted  by  Sathe  and  Tong  to  better  understand  the  heat 
and  fluid  flow  characteristics  inside  the  porous  medium.4 
At  the  present  time,  these  studies  use  an  assumed  value  for 
the  thermal  conductivity  of  the  Saran  Carbon  and  PCO  due  to 
lack  of  data.  To  generate  more  realistic  numerical  values 
the  assumed  values  should  be  replaced  by  experimentally 
obtained  values. 

Fourier’s  law  of  heat  conduction  was  used  in  obtaining 
the  thermal  conductivities  of  the  porous  material.  When 
applied  to  a  planar  medium  as  depicted  in  Figure  1, 

Fourier’s  law  yields  the  following  relationship 

k  =  QAx/AAT  (1) 

where  Q  is  the  energy  transferred  across  the  sample,  Ax  is 
the  sample  thickness,  A  is  the  area  perpendicular  to  the 
energy  transfer,  and  AT  is  the  temperature  difference 
between  the  hot  and  cold  surfaces  of  the  sample,  i.e.  Th  - 
Tc .  The  above  equation  is  valid  only  for  the  case  of  one¬ 
dimensional,  steady-state  energy  transfer,  and  a  non¬ 
temperature  dependent  (constant)  thermal  conductivity. 
Although  the  thermal  conductivity  of  materials  is  usually 
temperature  dependent,  the  dependency  effect  is  minimized  if 
AT  is  small.  The  thermal  conductivity  can  be  determined 
once  the  variables  on  the  right-hand  side  of  the  equation 


are  measured. 


Fig.  1  One-dinens ionel  heat  conduction  through  a  planar 
aediua 


CHAPTER  2 


EXPERIMENTAL  METHOD 

As  mentioned  ealier,  the  two  types  of  porous  material 
used  in  the  adsorption  compressors  of  interest  are  PCO  and 
Saran  Carbon.  To  determine  the  thermal  conductivities  of 
these  materials,  two  different  test  sections  were  designed 
and  built.  Two  test  sections  were  needed  because  of  the 
composition  of  the  available  porous  material.  The  PCO  was 
in  powder  form  while  the  Saran  Carbon  was  in  the  form  of  a 
rigid  rod.  Both  test  sections  produce  steady-state,  one¬ 
dimensional  heat  transfer  across  the  material  of  interest. 

2.1  Heat  Transfer  Considerations  for  the  Design  of  the  PCO 

Test  Apparatus 

To  measure  the  thermal  conductivity  of  the  PCO  a 
circular  guarded-hot-plate  apparatus5*6  was  designed  and 
built.  Figure  2  shows  the  basic  construction  of  a  guarded- 
hot-plate  apparatus.  The  bottom  heater  consists  of  two 
sections,  the  inner  or  test  section,  and  a  side  guard 
section.  To  minimize  heat  loss  in  the  lateral  direction 
from  the  test  section,  the  temperature  of  the  outer  hot 
plate  is  matched  with  that  of  the  inner  hot  plate  by 
appropriately  adjusting  the  side  guard  heater.  Similarly, 
the  rear  guard  heater,  the  top  plate,  and  the  insulation  are 
used  to  reduce  heat  transfer  in  the  upward  direction.  This 
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Cross-sectional  view  of  the  guarded-hot-plate 
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ensures  that  heat  transfer  from  the  test  section  can  be 
regarded  as  downward  and  one-dimensional  through  the  control 
volume  (C.V.).  To  arrive  at  the  dimensions  of  the 
apparatus,  a  heat  transfer  analysis  was  performed.  The 
purpose  was  to  determine  the  dimensions  of  the  different 
sections  so  that  heat  transfer  through  the  control  volume 
could  in  fact  be  treated  as  one-dimensional. 

A  Fortran  algorithm  SIMPLE  was  used  to  thermally  model 
the  idealized  guarded-hot-plate  apparatus,  as  shown  in 
Figure  3.  Adiabatic  condition  was  assumed  for  the  surface 
above  the  heaters.  The  heaters  were  modeled  as  Inconel 
(0.64  cm  thick)  discs  with  uniform  heat  generation.  The 
inner  and  outer  hot  plates  were  assumed  to  be  copper  (0.64 
cm  thick).  The  support  wall  (0.32  cm  thick)  and  the 
reference  plate  (0.64  cm  thick)  were  taken  as  stainless- 
steel.  SIMPLE  is  a  finite-difference  algorithm7  which  can 
be  used  to  solve  for  the  steady-state  temperature 
distribution  in  two-dimensional  objects.  Calculations  can 
be  performed  using  one  of  three  different  coordinate 
systems:  r-x,  x-y,  and  r-©.  The  coordinate  system  used  in 
this  study  was  r-x. 

The  temperature  distribution  calculated  in  the 
apparatus  can  be  used  to  provide  an  energy  balance  on  the 
control  volume  shown  in  Figure  3.  Since  the  amount  of 
energy  passing  through  the  sides  and  bottom  of  the  control 
volume  can  be  calculated  the  energy  loss  then  becomes  a 
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Fig.  3  Cross-sectional  view  of  the  idealized  guarded 
-hot-plate  apparatus 
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known  quantity.  The  energy  loss  on  a  percentage  basis  is 
defined  as  the  amount  of  energy  transferred  through  the 
sides  of  the  control  volume  divided  by  the  energy  generated 
in  the  inner  heater.  Thus,  it  is  an  indicator  of  how  good 
the  one-dimensional  heat  flow  assumption  is  for  the 
experiment.  An  acceptable  energy  loss  of  5  percent  or  less 
was  used  as  a  constraint  on  the  design  process. 

Figure  4  shows  the  relationship  between  the  percent 
energy  loss  and  the  ratio  of  D0  to  for  different  sample 
thicknesses.  In  this  case,  D0  is  the  outer  diameter  of  the 
side  guard  heater,  and  D4  is  the  outer  diameter  of  the  inner 
heater.  The  results  are  for  k  =  0.05  W/m*C,  which  should 
yield  a  conservative  estimate  of  the  energy  loss  since  it  is 
rather  small.  As  indicated  in  the  figure,  the  energy  loss 
on  a  percentage  basis  decreases  as  D0/Di  increases  for 
different  sample  thicknesses.  Also,  thicker  samples  result 
in  higher  energy  loss.  The  figure  was  used  to  determine  the 
dimensions  of  the  guarded-hot-plate  apparatus  for  the  PCO 
experiment.  This  was  accomplished  by  finding  the  region  on 
Figure  4  where  the  energy  loss  would  be  5  percent  or  less, 
and  then  picking  a  point  somewhere  in  the  region.  In  this 
experiment,  a  sample  thickness  of  0.64  cm  and  a  D0/Dj  ratio 
of  1.45  was  selected  to  meet  the  energy  loss  constraint. 
According  to  Figure  4,  the  percent  energy  loss  would  be  4.6 


percent . 
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Figure  5  shows  the  percent  energy  loss  as  a  function  of 
thermal  conductivity  for  the  dimensions  chosen  above.  As 
indicated  in  the  figure,  the  energy  loss  decreases 
exponentially  as  the  thermal  conductivity  increases. 
Therefore,  if  the  actual  thermal  conductivity  of  the  PCO  is 
greater  than  0.05  W/m*C,  the  percent  energy  loss  would  be 
less  than  4.6  percent. 

2.2  PCO  Experimental  Apparatus 

The  PCO  experimental  apparatus,  which  can  be  seen  in 
Figure  6,  is  suitable  for  thermal  conductivity  measurements 
on  materials  up  to  temperatures  of  700*C.  The  vacuum 
chamber  allows  for  measurements  under  vacuum  conditions  in 
order  to  eliminate  gaseous  heat  conduction.  The  dimensions 
of  the  apparatus  are  as  indicated  in  the  figure.  The  wall 
that  encloses  the  PCO  is  a  stainless  steel  tube  with  holes 
drilled  through  the  sides  to  facilate  evacuation  of  the  PCO. 
The  PCO  experimental  apparatus  was  placed  inside  a  vacuum 
chamber  by  bolting  the  reference  plate  to  a  flange  located 
at  the  bottom  of  the  chamber. 

The  hot  plate  sections  were  fabricated  from  0.64  cm 
thick  copper  plate  and  the  reference  plate  (cold  plate)  from 
1.27  cm  stainless-steel.  An  inner  heater  and  guard  heater 
were  installed  between  the  bottom  copper  plates  and 
insulation.  A  0.64  cm  thick  copper  plate  was  placed  on  top 
of  the  insulation  and  another  heater  (the  rear  guard  heater) 


PERCENT  ENERGY  LOSS 
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Fig.  5  Percent  energy  loss  against  thermal 
conductivity 
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was  placed  on  top  of  the  plate.  An  upper  piece  of 
insulation  was  placed  on  the  rear  guard  heater  and  a  0.64  cm 
thick  stainless-steel  plate  was  placed  on  top  of  the 
insulation. 

Eight  bare-wire  Chromel-Alumel  thermocouples,  0.038  cm 
diameter,  were  installed  in  grooves  machined  in  the  upward- 
face  of  the  copper  plates.  The  thermocouple  beads  were 
installed  in  small  holes  drilled  into  the  ends  of  the 
grooves.  The  bare  wires  were  coated  with  a  high  temperature 
cement  (Omega  CC)  and  placed  in  the  grooves.  This  ensured 
good  electrical  isolation  of  the  bare  wires  from  each  other 
as  well  as  from  the  copper  plates.  The  grooves  were  filled 
in  with  a  highly  conductive  material  (Steel  Seal)  to  ensure 
an  uniform  temperature  distribution  throughout  the  plates. 
Figure  7  shows  the  locations  of  the  thermocouples. 

Two  insulated  Chromel-Alumel  thermocouples  were 
installed  into  the  bottom  of  the  reference  plate  at  the  mid¬ 
plane  position.  The  thermocouples  were  held  into  place  by 
retaining  clips  attached  to  the  reference  plate.  The 
insulation  was  covered  with  stainless-steel  braiding  to 
protect  the  thermocouples  from  the  aluminum  oxide  sand  used 
in  the  constant  temperature  bed  which  will  be  dicussed 
later.  The  locations  of  the  reference  plate  thermocouples 
can  be  found  in  Figure  8. 

The  temperature  of  the  copper  plates  was  regulatad  by 
three  Digi-Sense  proportional  controllers  (model  #  2168-80). 


16 


The  temperature  controllers  automatically  controlled  the 
temperature  by  regulating  the  electrical  current  into  the 
heaters,  depending  upon  the  signal  they  received  from  the 
thermocouple  sensors  located  inside  the  copper  plates.  A 
power  pulse  meter  was  used  to  measure  the  amount  of  average 
power  delivered  to  the  inner  heater  at  steady-state 
conditions.  This  was  done  by  measuring  the  number  of  pulses 
delivered  to  the  inner  heater  in  a  certain  time  interval, 
and  the  resistance  as  well  as  the  r.m.s.  value  of  the 
voltage  drop  across  the  main  heater  (see  Appendix  A).  The 
amount  of  average  power  delivered  to  the  inner  heater  at 
steady-state  corresponds  to  Q  in  equation  (1). 

The  thermocouple  and  heater  leads  inside  the  chamber 
were  connected  to  a  vacuum  feedthrough  system  located  on  the 
side  of  the  chamber  (see  Figure  6).  The  feedthrough  was 
used  to  insure  no  leakage  into  the  chamber  from  the 
surrounding  room  air.  The  feedthrough  leads  were  connected 
to  a  thermocouple  switchbox  located  outside  the  chamber. 

The  error  introduced  by  the  feedthrough  system  was 
determined  by  calibrating  a  thermocouple  circuit  with  the 
feedthrough  system  installed  against  a  thermocouple  circuit 
without  the  feedthrough.  Data  taken  at  various  temperatures 
revealed  negligible  error  due  to  the  feedthrough  system  (see 
Appendix  B )  . 

Figure  9  shows  a  schematic  of  the  overall  system  used 
in  the  experiment.  Because  measurements  with  the  reference 
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plate  temperature  as  high  as  600*C  were  needed,  a  fluidized 
sand  bed  was  used  to  control  the  reference  plate  temperature 
by  immersing  the  vacuum  chamber  (with  the  reference  plate 
attached)  into  the  bed.  The  bed  was  capable  of  providing 
rapid  heat  transfer  and  precise  temperature  control  of  the 
reference  plate,  in  the  range  of  50*C  to  600*C.  The  bed 
heater  temperature  controller  was  used  to  maintain  the  bed 
temperature  to  within  ±  0.2*C  of  desired  setting.  A  FLUKE 
multimeter  (model  #  8842A)  was  used  to  collect  the  millivolt 
readings  from  the  thermocouples  installed  in  the 
experimental  apparatus. 

The  schematic  diagram  in  Figure  10  shows  the  major 
components  of  the  vacuum  system.  The  mechanical  roughing 
pump  was  used  to  pump  down  the  system  to  the  pressure  where 
the  turbomolecular  pump  could  start  pumping.  The 
turbomolecular  pump  was  used  to  obtain  the  desired  vacuum 
level.  The  thermocouple  and  ionization  gauges  were  used  to 
register  the  pressure  in  the  system  from  atmosphere  to  10~9 
torr.  A  variable  leak  valve  provided  a  means  to  control  the 

level  of  vacuum  inside  the  chamber.  For  the  PCO  powder  (10 

mjb  in  diameter)  a  vacuum  level  of  lxlO'3  torr  or  better  was 

needed.  This  insured  that  the  mean  free  path  of  the 

interstitial  fluids  would  exceed  the  spacing  between  the  PCO 
particles,  thereby  eliminating  gaseous  conduction.®  For 
most  of  the  experimental  runs,  a  vacuum  of  better  than  8x10“ 


5  torr  was  achieved. 


THERMOCOUPLE  GAUGE 
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Fig.  10  Conponents  of  the  vacuus  system 
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2.3  Experimental  Procedure 

PCO  was  placed  inside  the  apparatus  between  the  bottom 
copper  plates  and  the  reference  plate  (see  Figure  6).  This 
was  done  by  first  attaching  the  stainless-steel  tube  to  the 
reference  plate.  Then,  by  using  a  ram  assembly,  loose  PCO 
powder  (approximately  10  pm  in  diameter)  was  compressed,  by 
JPL,  to  a  depth  of  0.64  cm  and  density  of  3.5  grams/cm3, 
into  the  bottom  of  the  stainless-steel  tube.  The  hot  plate 
section  was  then  placed  on  top  of  the  PCO  and  the  reference 
plate  was  bolted,  after  attaching  the  thermocouple  and 
heater  leads  to  the  feedthrough  system,  to  the  vacuum 
chamber  as  explained  in  section  2.2. 

Once  the  PCO  apparatus  had  been  bolted  to  the  vacuum 
chamber  the  system  was  placed  into  the  constant  temperature 
fluidized  bed.  The  vacuum  system  was  turned  on  and  the 
vacuum  level  inside  the  chamber  reached  8xl0~5  torr.  After 
outgassing,  all  temperature  controllers  were  set  to  the 
desired  AT,  typically  around  25*C,  across  the  top  and  bottom 
surfaces  of  the  PCO.  The  level  of  vacuum  inside  the  chamber 
was  then  set  to  a  predetermined  reading  and  data  were  taken 
at  30  minute  intervals  until  steady-state  conditions  were 
obtained.  In  this  experiment,  steady-state  conditions 
existed  when  the  desired  AT  did  not  vary  more  than  t  0.5*C 
over  a  60-min  interval.  The  procedure  was  then  repeated  for 
different  hot-  and  reference-  plate  temperatures.  The  data 
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collected  consisted  of  the  temperature  at  each  thermocouple 
location,  number  of  pulses  to  the  inner  heater  in  a  time 
interval,  and  bed  temperature.  The  thermal  conductivity  was 
then  calculated  using  equation  (1).  The  uncertainties  in 
the  values  of  the  thermal  conductivity  ranged  from  t  2.44 
percent  to  t  2.56  percent  as  given  in  Appendix  C. 

2.4  Heat  Transfer  Considerations  for  the  Design  of  the 

Carbon  Test  Section 

To  measure  the  thermal  conductivity  of  the  Saran  Carbon 
a  different  type  of  apparatus,  based  on  previous 
designs9*10'11,  was  designed  and  built.  The  basic  apparatus 
can  be  seen  in  Figure  11.  The  apparatus  consists  of  two 
columns  with  the  test  specimen  sandwiched  between  them.  The 
objective  is  to  achieve  one-dimensional  heat  conduction  from 
one  end  of  the  carbon  rod  to  the  other.  A  heat  source  is 
attached  to  the  top  column  to  provide  the  needed  heat 
transfer  into  the  test  specimen. 

SIMPLE  was  used  to  thermally  model  the  apparatus.  This 
was  done  to  determine  if  heat  transfer  in  the  carbon  rod 
could  be  treated  as  one-dimensional.  The  same  procedure  as 
discussed  before  (see  section  2.1)  was  used  to  determine  the 
energy  loss  through  the  sides  of  the  control  volume  shown  in 
Figure  11.  In  this  investigation,  these  losses  never 
exceeded  1.6  percent,  indicating  the  one-dimensional  heat 
conduction  assumption  was  acceptable. 


HEAT  SOURCE 


C.V. 


BOTTOM 

COLUMN 


Fif.  11  Basic  apparatus  used  for  the  Saran  Carbon 
exper iaent 
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2.5  Saran  Carbon  Experimental  Apparatus 

The  Saran  Carbon  experimental  apparatus,  which  can  be 
seen  in  Figure  12,  is  suitable  for  thermal  conductivity 
measurements  on  rigid  rod  samples  in  the  temperature  range 
between  -20*C  and  200*C.  The  dimensions  of  the  apparatus 
are  as  indicated  in  the  figure.  The  apparatus  was  composed 
of  a  top  aluminum  column  with  a  heater  attached,  a  spacer 
made  of  insulation,  a  load  plate  with  a  support  structure,  a 
bottom  stainless-steel  column,  and  a  stainless  steel 
reference  plate.  The  carbon  was  sandwiched  between  the  top 
and  bottom  columns.  The  load  plate  was  used  to  insure  good 
surface  contact  between  the  columns  and  the  carbon  by 
providing  an  axial  load  onto  the  column  as  the  plate  was 
bolted  to  the  support  structure.  The  insulation  was  used  to 
minimize  heat  transfer  into  the  load  plate.  A  240-W,  120 
Volt  coil  heater  was  wrapped  around  the  top  aluminum  column 
to  provide  the  needed  heat  source.  A  highly  reflective 
radiation  shield  was  installed  around  the  column  to  reduce 
radial  heat  losses  due  to  radiation  heat  transfer. 

Ten  insulated  Chromel-Alumel  thermocouples,  0.099  cm 
diameter,  were  installed  in  holes  drilled  into  the 
centerline  position.  All  of  the  thermocouples  were  held 
into  place  by  a  highly  conductive  epoxy  (OMBGABOND  200). 
Using  the  known  values  of  the  temperatures  and  distances 
between  the  thermocouple  locations  in  the  stainless  steel 


THERMOCOUPLE  LOCATIONS 
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Fig.  12  Saran  Carbon  experimental  apparatus 
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column,  as  well  as  the  thermal  conductivity  of  stainless 
steel,  the  energy  transfer,  Q,  into  the  stainless  steel 
column  could  be  determined.  With  Q  as  a  known  quantity 
equation  (1)  could  be  used  to  solve  for  the  thermal 
conductivity  of  the  Saran  Carbon. 

The  temperature  of  the  apparatus  was  controlled  by  one 
Digi-Sense  proportional  temperature  controller  as  discussed 
in  section  2.2.  One  of  the  installed  thermocouples  provided 
the  needed  temperature  sensor  for  the  controller.  The 
temperature  of  the  reference  plate  was  controlled  by  either 
the  fluidized  bed  or  a  constant  temperature  bath,  depending 
on  the  desired  temperature  of  the  plate.  The  constant 
temperature  fluidized  bed  provided  control  between  50*C  and 
200*C,  and  the  constant  temperature  bath  provided  control 
between  -20*C  and  49*C. 

The  thermal  conductivity  of  the  carbon  was  determined 
in  a  vacuum  environment.  This  was  necessary  to  eliminate 
conduction  losses  to  the  surrounding  air  and  also  to 
eliminate  the  effects  of  any  interstitial  fluids  on  the 
thermal  conductivity.  For  the  Saran  Carbon  powder  (2  /im  in 
diameter)  a  vacuum  level  of  lxlO-3  torr  or  better  was 
necessary  to  eliminate  gaseous  conduction  as  discussed 
previously  in  section  2.2.  The  reference  plate  was  bolted 
to  the  vacuum  chamber.  The  thermocouple  and  heater  leads 
were  passed  from  the  inside  of  the  chamber  to  an  outside 
switchbox  as  discussed  previously  in  section  2.2. 
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2.6  Experimental  Procedure 

The  Saran  Carbon  test  specimen  was  formed,  by  JPL,  into 
an  approximately  0.127  cm  OD  rod  by  compressing  the  carbon 
powder  into  a  die  assembly  to  a  density  of  1.135  grams/cm3, 
and  a  porosity  of  47.2  percent.  After  compression,  the  rod 
was  removed  from  the  die  assembly  and  cut  to  a  length  of 
2.06  cm.  The  test  column  was  fully  assembled  and  bolted, 
after  attaching  the  inner  leads  to  the  vacuum  feedthrough 
system  and  installing  the  radiation  shield,  to  the  vacuum 
chamber . 

After  assenbling  the  apparatus,  the  vacuum  chamber  was 
immersed  into  the  fluidized  bed  or  the  constant  temperature 
bath  as  mentioned  before.  The  vacuum  system  was  turned  on 
and  allowed  to  reach  a  constant  level  of  vacuum  inside  the 
chamber.  For  all  of  the  experimental  runs  a  vacuum  of 
better  than  lxlO-4  torr  was  achieved.  After  outgassing,  all 
temperature  controllers  were  set  to  the  desired  AT  across 
the  carbon  specimen.  Data  was  then  taken  at  30  minute 
intervals  until  steady-state  conditions  were  insured  in  the 
same  manner  as  discussed  in  section  2.3.  The  data  consisted 
of  the  temperatures  at  each  thermocouple  location,  vacuum 
level  inside  the  chamber,  bath  temperature,  and  power 
setting  of  the  heater.  The  AT  across  the  carbon  was  then 
adjusted  to  a  different  setting  and  the  procedure  was 
repeated  once  again.  The  uncertainties  in  the  values  of  the 
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thermal  conductivity  ranged  from  t  5.43  percent  to  t  12.94 
percent  as  given  in  Appendix  D.  For  approximately  eighty 
percent  of  the  runs  the  uncertainty  in  the  values  of  the 
thermal  conductivity  was  less  than  8.17  percent. 


CHAPTER  3 


RESULTS  AND  DISCUSSION 


Table  1  shows  a  summary  of  the  test  conditions  and 
results  for  the  PCO  experiment.  In  this  case,  the  mean 
temperature  (Tmean)  is  the  arithmetic  nean  between  the  hot 
and  cold  plate  temperatures.  The  thermal  conductivity 
values  were  calculated  from  equation  (1)  as  mentioned  in 
section  2.2. 


Table  1  PCO  experimental  conditions  and  results 


Run 

* 

T 

mean  * 

•c 

Th* 

*C 

T 

» 

*C 

Pulse , 
count 

At, 

secs 

k, 

W/m*  C 

1 

312.57 

326. 15 

299. 00 

7005.0 

5212.40 

0.1471 

2 

323.40 

335.20 

311.60 

5830.0 

5340.32 

0.1378 

3 

362.00 

375.50 

348.50 

9274.0 

5232.58 

0.1957 

4 

410.15 

422.90 

397.40 

12681.0 

5099.58 

0.2867 

5 

457.60 

470.30 

444.90 

18200.0 

4907.68 

0.4207 

6 

507.65 

519.65 

495.65 

23020.0 

5335.63 

0.5211 

7 

556.43 

568.20 

544.65 

27990.0 

5394.38 

0.6389 

8 

591.95 

614.60 

569.30 

49466.0 

4080.28 

0.7624 

Figure  13  illustrates  the  results  of  the  PCO 
experiment.  As  shown,  the  thermal  conductivity  of  the  PCO 
has  a  strong  temperature  dependency.  The  line  in  the  figure 
was  obtained  from  the  following  equation 

k  =  -0.00013  -0.319659E-3TBean  +  0 . 160629E-8TBe a „3  (2) 

which  was  a  least-squares-error  fit  of  the  data  with  TBean 
in  Kelvin.  The  error  sum  of  squares  is  0.14  percent, 


indicating  a  very  good  fit  through  the  points. 


The  first 
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two  terms  on  the  right-hand  side  allows  for  a  linear 
temperature  variation.  The  cubic  term  is  included  to 
account  for  possible  radiation  contribution  to  the  thermal 
conductivity.  The  cubic  term  arises  from  the  basic 
definition  of  the  total  thermal  conductivity  of  a  porous 
material.8  In  this  case,  it’s  magnitude  is  comparable  to 
the  other  terms  in  the  equation,  indicating  an  appreciable 
radiation  contribution  to  the  thermal  conductivity. 

Equation  (2)  was  used  to  calculate  the  thermal 
conductivity  at  each  mean  temperature,  and  then  these  values 
were  compared  with  the  experimental  results  as  shown  in 
Table  2.  As  indicated  in  the  table,  the  greatest  difference 
between  the  correlation  and  the  experimentally  obtained 
values  is  8.68  percent,  with  75  percent  of  the  results 
within  5.84  percent  or  less. 


Table  2  Comparison  of  calculated  values  from  equation  (2) 
with  experimentally  obtained  values  of  the  PCO 


T 

mean  * 

*  c 

Calculated  k, 
W/m*  C 

Experimental  k, 
W/m*  C 

Difference, 

percent 

312.57 

0. 1352 

0.1471 

8.68 

323.40 

0. 1499 

0.1378 

8.65 

362.00 

0.2081 

0.1957 

5.84 

410. 15 

0.2935 

0.2867 

2.22 

457.60 

0.3927 

0.4207 

6.94 

507.65 

0.5145 

0.5211 

1.26 

556.43 

0.6512 

0.6389 

1 . 87 

591.95 

0.7626 

0.7624 

0.85 

Table  3  gives  a  summary  of  the  test  conditions  and 
results  for  the  Saran  Carbon  experiment.  In  this 


experiment,  the  mean  temperature  is  the  arithmetic  mean 
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between  thermocouples  #  4  and  #  6  located  in  the  carbon  rod 
(see  Appendix  F).  Equation  (1)  was  used  to  obtain  the 
thermal  conductivity  values  shown  in  the  table.  This  was 
accomplished  in  the  following  manner.  First,  a  plot  of 
temperature  against  thermocouple  location  was  generated  for 
the  stainless  steel  column.  Then,  a  straight  temperature 
distribution  line  was  fit  through  the  three  data  points  and 
was  used  to  determine  Q  through  the  stainless  steel  column. 
To  satisfy  energy  balance,  heat  transfer  through  the  carbon 
rod  must  have  the  same  Q.  Similarly,  a  straight  line  was 
fit  through  the  three  temperature  data  points  for  the  carbon 
rod.  The  thermal  conductivity  of  the  carbon  was  solved  for 
by  equating  Q/kBEA  to  the  slope  of  the  straight  line. 


Table  3  Saran  Carbon  experimental  conditions  and  results 


Run , 

# 

T 

Am  e  a  n  * 

'C 

Th. 

'C 

Tc* 

•c 

k, 

W/m*  C 

1 

2.89 

8.20 

-2.40 

0.785 

2 

33.54 

46.90 

20. 17 

0.702 

3 

66.22 

87.81 

44.63 

0.709 

4 

77.58 

85.23 

69.93 

0.883 

5 

78.64 

86.28 

70.98 

0.878 

6 

115.77 

130.12 

101.42 

0.976 

7 

133.00 

150.41 

115.60 

1.040 

8 

133.10 

152.21 

114.00 

1.077 

9 

133.83 

141.96 

125.70 

1.082 

10 

142.16 

160.72 

123.60 

1.100 

11 

142.80 

162.37 

123.23 

1.106 

12 

166.72 

189.99 

143.44 

1.142 

13 

193.70 

211 . 96 

175.42 

1.166 

F i gure 

14  shows  the 

experimental 

thermal 

conductivity 

plotted  against  mean  temperature.  The  data  was  obtained 
over  the  temperature  range  used  in  the  Saran  Carbon 


33 


adsorption  compressors.  The  solid  symbols  represent  the 
data  obtained  using  the  sand  bed  and  the  open  symbols 
represent  the  data  obtained  using  the  constant  temperature 
bath.  As  shown  in  the  figure,  the  thermal  conductivity 
increases  with  a  corresponding  increase  in  the  mean 
temperature.  The  following  equation 

k  =  0.20359  +  0. 154475E-2TBean  +  0 . 303628E-9TBe a n 3  (3) 

was  again  obtained  through  a  least-squares-error  fit  of  the 
data  with  Tmean  in  Kelvin.  For  this  fit,  the  error  sum  of 
squares  is  4.28  percent.  In  this  case,  the  contribution  of 
radiation  is  negligible  since  the  cubic  term  is  very  small 
compared  to  the  first  two  terms,  even  for  the  data  point 
with  the  highest  mean  temperature.  In  fact,  a  fit  with  a 
linear  function  yielded 

k  =  -0.10115  +  0.280980E-2TBean  (4) 

with  4.41  percent  as  the  error  sum  of  squares  and  TBean  in 
Kelvin.  The  line  in  Figure  14  was  obtained  from  equation 
(3).  It  is  clear  that  a  linear  variation  is  sufficient  to 
allow  for  the  temperature  dependency. 

Equation  (4)  was  used  to  calculate  the  thermal 
conductivity  at  each  mean  temperature  value,  and  then  these 
values  were  compared  with  the  experimental  results  as  shown 
in  Table  4.  As  indicated  in  the  table,  the  greatest 
difference  between  the  calculated  values  and  the 
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experimentally  obtained  values  is  16.78  percent,  with  84.6 
percent  of  the  results  within  7.64  percent  or  less. 


Table  4  Comparison  of  calculated  values  from  equation  (4) 
with  experimentally  obtained  values  of  the  Saran 
Carbon 


T 

mean  * 

*c 

Calculated  k, 
W/ra*  C 

Experimental  k, 
W/m*  C 

Difference, 

percent 

2.89 

0.6740 

0.7850 

16.47 

33.54 

0.7601 

0.7020 

7.64 

66.22 

0.8520 

0.7090 

16.78 

77.58 

0.8840 

0.8830 

0.11 

78.64 

0.8870 

0.8780 

1.01 

115.77 

0.9912 

0.9760 

1.54 

133.00 

1 . 0396 

1.0400 

0.04 

133.10 

1.0399 

1.0770 

3.56 

133.83 

1.0419 

1 . 0820 

3.85 

142. 16 

1.0654 

1.1000 

3.24 

142.80 

1.0672 

1.1060 

3.63 

166. 72 

1.1344 

1.1420 

0.67 

193.70 

1.2102 

1.1660 

3.65 

The  thermal  conductivity  of  the  Saran  Carbon  is  further 
compared  to  the  thermal  conductivity  of  porous  carbon,  which 
is  assumed  to  be  given  by  the  following  relationship 

k  =  +  ( 1-e ) ks  (5) 

where  e  is  the  porosity  of  the  porous  carbon,  which  is 
assumed  to  be  the  same  as  the  Saran  Carbon,  kt  is  the 
thermal  conductivity  of  the  interstitial  fluid  (in  this 
experiment  =0),  and  ks  is  the  thermal  conductivity  of 
solid  carbon.  Table  5  shows  the  comparison  between  the 
Saran  Carbon  and  porous  carbon  at  different  mean  temperature 
values  with  t  -  47.2  percent.  The  largest  difference  is 
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24.70  percent,  with  84.6  percent  of  the  exper imental  results 
falling  within  9.88  percent  of  each  other. 


Table  5  Comparison  of  Saran  Carbon  with  solid  carbon 


T 

mean  * 

*C 

Th, 

•c 

Tc. 

•  C 

k, 

W/m*  C 

^B  » 

W/m*  C 

(l-i)k,  , 

W/m*  C 

2.89 

8.20 

-2.40 

0.785 

1.449 

0.765 

33.54 

46.90 

20.17 

0.702 

1.558 

0.822 

66.22 

87.81 

44.63 

0.709 

1.674 

0.883 

77.58 

85.23 

69.93 

0.883 

1.715 

0.905 

78.64 

86.28 

70.98 

0.878 

1.718 

0.907 

115.77 

130.12 

101.42 

0.976 

1.850 

0.976 

133.00 

150.41 

115.60 

1.040 

1.900 

1.002 

133.10 

152.21 

114.00 

1.077 

1.901 

1.002 

133.83 

141.96 

125.70 

1.082 

1.902 

1.003 

142. 16 

160.72 

123.60 

1.100 

1.913 

1.010 

142.80 

162.37 

123.23 

1.106 

1.915 

1.010 

166.72 

189.99 

143.44 

1.142 

1.950 

1 . 029 

193.70 

211.96 

175.42 

1.166 

1.990 

1.050 

CHAPTER  4 


CONCLUSIONS 


Numerical  studies  are  being  conducted  by  Sathe  and  Tong 
to  better  understand  the  heat  and  fluid  flow  characterics 
inside  porous  adsorbents4  used  in  laboratory  adsorption 
compressors.  These  compressors  use  Praseodymium-Cerium- 
Oxide  (PCO)  and  Saran  Carbon  as  the  adsorbents.  At  the 
present  time,  these  studies  use  an  assumed  value  for  the 
thermal  conductivity  of  the  PCO  and  Saran  Carbon  due  to  lack 
of  data.  It  was  the  objective  of  the  present  investigation 
to  experimentally  obtain  values  for  the  thermal  conductivity 
of  both  materials. 

Since  the  porous  materials  came  in  different  forms,  the 
PCO  in  a  powder  and  the  Saran  Carbon  in  a  rigid  rod,  two 
types  of  test  sections  were  needed  for  the  thermal 
conductivity  measurements.  A  circular  guarded-hot-plate 
apparatus  was  designed  and  built  for  the  PCO  experiment, 
while  a  column  apparatus  was  used  for  the  Saran  Carbon. 

Both  test  sections  simulate  steady-state,  one-dimensional 
heat  transfer  across  the  material  of  interest. 

Results  of  the  thermal  conductivity  measurements 
suggest  a  strong  temperature  dependency  for  both  materials. 

A  least-squared-error  fit  program  was  used  to  fit  a  line 
through  the  PCO  and  Saran  Carbon  data  points.  Comparison 
between  the  calculated  values  from  equation  (2)  and  the 
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experimentally  obtained  values  of  the  PCO  thermal 
conductivity  indicates  that  75  percent  of  the  data  is  within 
5.84  percent  or  less  of  the  line  in  Figure  13. 

Comparison  between  the  Saran  Carbon  and  porous  carbon  shows 
the  greatest  difference  as  24.7  percent,  with  84.6  percent 
of  the  experimental  results  within  9.88  percent  of  the 
porous  values  (see  Table  5). 

Radiation  effects  significantly  contributed  to  the 
thermal  conductivity  measurements  of  the  PCO,  while  showing 
a  negligible  contribution  to  the  thermal  conductivity  of  the 


Saran  Carbon. 
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Power  Pulse  Meter 
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The  power  pulse  meter  is  a  device  that  counts  the 
number  of  pulses  delivered  to  the  inner  heater  in  a  certain 
time  interval.  The  meter  senses  the  half  wavelength  of  an 
a.c.  power  pulse.  When  full  power  is  applied  the  meter 
counts  120  pulses  per  second,  because  there  are  120  half 
wavelengths  in  a  60  Hz  power  supply. 

For  example,  consider  the  following  case.  The 
resistance  (R)  of  the  inner  heater  is  116.58  Q.  At  full 
power  the  r.m.s.  voltage  drop  (Ve)  across  the  inner  heater 
was  117.05V.  Thus, 

Full  power  =  Ve2/R  =  117. 52W. 

Now,  with  2448  power  pulses  supplied  to  the  inner 
heater  in  a  time  interval  (At)  of  392.72  secs  we  have: 

Average  power  =  Full  power  x  2448/(120  x  392.72)  =  6.1W 
Therefore,  average  power  corresponds  to  Q  in  equation  (1). 


APPENDIX  B 

Feedthrough  Calibration 
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In  both  experiments  a  feedthrough  system  was  used  to 
collect  mv  reading  from  various  thermocouples  located  inside 


the  vacuum  chamber.  To  verify  if  any  error  would  be 
introduced  by  the  feedthrough  system  two  thermocouple 
circuits  for  each  experiment  were  built.  One  circuit  had 
the  feedthrough  system  installed  (WFT)  while  the  other  did 
not.  Data  was  taken  for  each  circuit  across  a  range  of 
temperatures  (see  Table  6)  and  then  the  circuit  with  the 
feedthrough  was  calibrated  against  the  circuit  without  the 
feedthrough . 


Table  6  Calibration  of  thermocouples 


Temp , 

•c 

Saran 

WFT, 

mv 

Carbon 

Other , 

mv 

PCO 

WFT, 

mv 

Other , 
mv 

80 

3.274 

3.268 

3.283 

3.277 

100 

4.099 

4.096 

4.105 

4.103 

120 

4.923 

4.921 

4.931 

4.  929 

140 

5.740 

5.744 

5.749 

5.746 

160 

6.586 

6.581 

6.572 

6.574 

200 

8.230 

8.232 

8.219 

8.217 

230 

9.468 

9.462 

9.425 

9.418 

270 

11.096 

11.089 

11.078 

11.086 

310 

12.753 

12.748 

12.718 

12.723 

350 

14.385 

14.380 

14.350 

14.354 

400 

16.426 

16.423 

16.398 

16.399 

450 

18.463 

18.460 

18.438 

18.440 

500 

20.494 

20.491 

20.482 

20.485 

APPENDIX  C 

Error  Analysis  for  the  PCO 
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From  equation  (1)  we  have: 
k  =  QAx/AAT 

but  Q  =  Ve2 ( count ) /120RAt 

where  count  is  the  number  of  pulses  delivered  to  the  inner 
heater  in  a  time  interval  (At),  and  Ve2  ,  R  are  defined  in 
Appendix  A.  So  k  is  now: 

k  =  ( Ve2 count Ax/ 1 20RAt ) /AAT 
The  root  mean  square  error  is  given  by: 

Ak/k  =  ^(2AVe/AVe)2  +  ( Acount/count ) 2  +  (AR/R)2  +(A(At)/At)2 
+  ( A ( Ax ) / Ax ) 2  +  (AA/A)2  +  (ATh/AT)2  +  ( ATC /AT ) 2 f 1 / 2 
For  example  consider  the  following  case: 


AVe 

= 

1.57V 

Ve 

= 

117.05V 

Acount 

= 

0.5 

count 

= 

2448 

AR 

= 

0.005  Q 

R 

= 

116.58  0 

A  (  At  ) 

= 

0.005  secs 

At 

= 

392.72  secs 

A  (Ax) 

= 

1.016E-4  m 

Ax 

- 

6 . 35E-3  m 

AA 

= 

1.805E-6  m 

A 

= 

6.43E-3  m 

*Th 

= 

0.2  *C 

AT 

= 

25  *C 

ATC 

= 

• 

CM 

O 

The  error  in  the  resistance,  and  voltage  is  from  the 
FLUKE  8842A  multimeter  used  in  the  experiment. 

Substituting  the  above  values  in  the  root  mean  square 
error  equation  yields: 

Ak/k  =  0.05044  i.e.  an  uncertainty  of  t  5.044  percent. 


APPENDIX  D 

Error  Analysis  for  the  Saran  Carbon 
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Again,  from  equation  (1)  we  have: 
k  =  QAxcar/AATcar 
but  for  this  case, 

Q  =  kBBAATBB/AxBB 

where  ks s  is  the  thermal  conductivity  of  the  stainless  steel 
column,  ATS s  is  the  temperature  difference  ( T3  -  Tx )  in  the 
stainless  steel  column,  Axss  is  the  distance  between 
thermocouple  locations  in  the  stainless  steel  column,  and 
ATcar  is  the  temperature  difference  (T6  -  T4 )  in  the  Saran 
Carbon-  So  k  is  now: 

k  =  (kasAATBB/Ax8B  )Axcar/AATcar 
The  root  mean  square  error  is  given  by: 

Ak/k  =  \  (AkBB/kBS  )2  +  (  2  ( AA/ A) 2  )  +(ATh/ATBB)2  +(ATC/ATBB)2  + 
(A(AxBB)/AxBB)2  +  (A(Ax)/Axcar)2  +  (ATh/ATcar)2  +  (ATc/ATcar)2 
Consider  the  following  example: 


=  0.800  W/m*  C 

kB» 

=  16.013  W/m*  C 

AA 

=  2.49E-7  m 

A 

=  1.224E-4  m 

j= 

H 

< 

=  ATC  =  0.2  *C 

AT 

‘si 

=  17.027  *C 

) 

=  1.27E-6  m 

Axss 

=  4.048E-2  m 

A(  Ax  ) 

=  1.27E-6  m 

i,(csr 

=  9.855E-3  m 

AT 

‘car 

=  38.145  *C 

where  Axc a r  is  the  distance  between  thermocouples  in  the 
Saran  Carbon. 


Substituting  the  above  values  in  the  root  mean  square 
error  equation  yields: 

Ak/k  =  0.05328  i.e.  an  uncertainty  of  5.328  percent. 


APPENDIX  E 

PCO  Experimental  Data 
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Table  7  PCO  experimental  data  for  run  #  1 


Thermocouple, 

# 

e  .  m  •  f* «  i 

nv 

*C 

1 

13.299 

326. 103 

3 

13.306 

326.296 

4 

12.450 

305.580 

6 

13.208 

323.908 

7 

13.230 

324.439 

9 

12.180 

299.033 

10 

12.173 

298.864 

Thermocouple  numbering  code  as  per  location: 

1-3  -  inner  copper  plate 

4-5  -  outer  copper  plate 

6-8  -  top  copper  plate 

9-10  —  reference  plate 

note  :  thermocouples  #  2,5,8  were  used  to  regulate 
controllers 


inner  heater  controller  temperature 
guard  heater  controller  temperature 
top  heater  controller  temperature 
bed  controller  temperature 
resistance  of  inner  heater 
count 
At 

vacuum 

r.m.s.  voltage  drop 
across  inner  heater 
Ax 

area  perpendicular  to 
heat  transfer 


325.2  *C 

307.2  *C 
337.1  *C 
305.0  *C 
84.20  Q 
7005 

5212.41  secs 
9E-5  torr 

122.35  volts 
6.35E-3  m 

3. 167E-3  m2 


Table  8  PCO  <  perimental  data  for  run  #  3 


Thermocouple , 

* 

C  •  D  •  f  •  i 
mv 

•c 

1 

15.332 

374.868 

3 

15.348 

375.250 

4 

14.380 

352.099 

6 

15.235 

372.553 

7 

15.268 

373.341 

9 

14.244 

348.837 

10 

14.218 

348.213 
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inner  heater  controller  temperature 
guard  heater  controller  temperature 
top  heater  controller  temperatue 
bed  controller  temperature 
resistance  of  inner  heater 
count 
At 

vacuum 

r.m.s.  voltage  drop 
across  inner  heater 


376.0  *C 

356.3  *C 

388.3  *C 

355.0  *C 
85.47  0 

9274 

5232.58  secs 
1.5E-4  torr 


122.63  volts 


Table  9  PCO  experimental  data  for  run  #  4 


Thermocouple , 

# 

e  .  m .  f .  , 
mv 

♦c 

1 

17.343 

422.615 

3 

17.365 

422.135 

4 

16.280 

397.432 

6 

17.318 

422.024 

7 

17.321 

422.095 

9 

16.292 

397.717 

10 

16.268 

397.147 

inner  heater  controller  temperature 
guard  heater  controller  temperature 
top  heater  controller  temperature 
bed  controller  temperature 
resistance  of  inner  heater 
count 
At 

vacuum 

r.m.s.  voltage  drop 
across  inner  heater 


425.6  *C 
401.5  *C 
441.4  *C 
405.0  *C 
85.47  0 
12681 

5099.58  secs 
1.4E-4  torr 


121.92  volts 


APPENDIX  F 

Saran  Carbon  Experimental  Data 
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Table  10  Saran  Carbon  experimental 

data  for  run  #  1 

Thermocouple,  e.m.f., 

« 

mv 

•c 

1 

-0.644 

-15.300 

2 

-0.590 

-14.000 

3 

-0.543 

-12.869 

4 

-0. 109 

-2.405 

5 

0. 108 

2.836 

6 

G.330 

8.204 

7 

0.901 

22.043 

8 

0.916 

22.407 

9 

0.920 

22.505 

Thermocouple  numbering  code  as  per  location: 

1  - 

3  -  bottom  stainless  steel  column 

4  - 

6  -  Saran  Carbon  sample 

7  - 

9  -  top  aluminum  column 

note  : 

thermocouple  #  10  was  used  to  regulate  controller. 
Thermophysical  properties  of  stainless  steel  were 

evaluated  using  the  temperature 
location  #  2 

at  thermocouple 

distance  between  thermocouple  #  1  and  # 

3  =  4.0481E-2  m 

distance  between  thermocouple  #  4  and  # 

6  =  9 . 8552E-3  m 

area  perpendicular  to  heat  transfer 

=  1.2240E-4  m2 

variac 

setting 

=  5 

constant  temperature  bath 

=  -22.0  *C 

controller  temperature 

=  22.5  *C 

Table  11  Saran  Carbon 

experimental 

data  for  run  #  4 

Thermocouple , 

6  •  D  •  f  •  f 

# 

mv 

•c 

1 

1.030 

25. 176 

2 

1.111 

27.144 

3 

1.182 

28.870 

4 

2.863 

69.928 

5 

3.150 

76.971 

6 

3.486 

85.227 

7 

6.406 

157.281 

8 

6.338 

155.600 

9 

6.442 

158. 195 

variac  setting 


16 
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constant  temperature  bath 
controller  temperature 


16.0  ‘C 
158.4  *C 
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